Arbitrary-Order Trigonometric Fourier Collocation Methods for Multi-Frequency Oscillatory Systems
نویسندگان
چکیده
We rigorously study a novel type of trigonometric Fourier collocation methods for solving multi-frequency oscillatory second-order ordinary differential equations (ODEs) q′′(t) + Mq(t) = f(q(t)) with a principal frequency matrix M ∈ Rd×d. If M is symmetric and positive semi-definite and f(q) = −∇U(q) for a smooth function U(q), then this is a multi-frequency oscillatory Hamiltonian system with the Hamiltonian H(q, p) = p p/2 + q Mq/2 + U(q), where p = q′. The solution of this system is a nonlinear multi-frequency oscillator. The new trigonometric Fourier collocation method takes advantage of the special structure brought by the linear term Mq and its construction incorporates the idea of collocation methods, the variation-of-constants formula and the local Fourier expansion of the system. The properties of the new methods are analysed. The analysis in the paper demonstrates an important feature, namely that the trigonometric Fourier collocation methods can be of an arbitrary order and when M → 0, each trigonometric Fourier collocation method creates a particular Runge–Kutta–Nyström-type Fourier collocation method, which is symplectic under some conditions. This allows us to obtain arbitrary high order symplectic methods to deal with a special and important class of systems of second-order ODEs in an efficient way. The results of numerical experiments are quite promising and show that the trigonometric Fourier collocation methods are significantly more efficient in comparison with alternative approaches that have previously appeared in the literature.
منابع مشابه
Analysis of a collocation method for integrating rapidly oscillatory functions
A collocation method for approximating integrals of rapidly oscillatory functions is analyzed. The method is efficient for integrals involving Bessel functions J,.(rx) with a large oscillation frequency parameter r, as well as for many other oneand multi-dimensional integrals of functions with rapid irregular oscillations. The analysis provides a convergence rate and it shows that the relative ...
متن کاملNumerical Integrators for Highly Oscillatory Hamiltonian Systems: A Review
Numerical methods for oscillatory, multi-scale Hamiltonian systems are reviewed. The construction principles are described, and the algorithmic and analytical distinction between problems with nearly constant high frequencies and with timeor state-dependent frequencies is emphasized. Trigonometric integrators for the first case and adiabatic integrators for the second case are discussed in more...
متن کاملA Family of Trigonometrically Fitted Enright Second Derivative Methods for Stiff and Oscillatory Initial Value Problems
A family of Enright’s second derivative formulas with trigonometric basis functions is derived using multistep collocation method. The continuous schemes obtained are used to generate complementary methods. The stability properties of the methods are discussed.Themethodswhich can be applied in predictor-corrector formare implemented in block formas simultaneous numerical integrators over nonove...
متن کاملNumerical solution of second-order hyperbolic telegraph equation via new cubic trigonometric B-splines approach
This paper presents a new approach and methodology to solve the second-order one-dimensional hyperbolic telegraph equation with Dirichlet and Neumann boundary conditions using the cubic trigonometric B-spline collocation method. The usual finite difference scheme is used to discretize the time derivative. The cubic trigonometric B-spline basis functions are utilized as an interpolating function...
متن کاملTrigonometric Cubic B-spline Collocation Method for Solitons of the Klein-Gordon Equation
In the present study, we derive a new B-spline technique namely trigonometric B-spline collocation algorithm to solve some initial boundary value problems for the nonlinear Klein-Gordon equation. In order to carry out the time integration with Crank-Nicolson implicit method, the order of the equation is reduced to give a coupled system of nonlinear partial differential equations. The collocatio...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید
ثبت ناماگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید
ورودعنوان ژورنال:
- Foundations of Computational Mathematics
دوره 16 شماره
صفحات -
تاریخ انتشار 2016